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ABSTRACT 

The  maximum  effect  of  random  errors  in  a five  station  firing  range 
is  investigated.  Because  the  data  reduction  method  for  drag  coefficient, 
K«,  requires  fitting  the  time-distance  data  with  a cubic  polynomial  by 

the  method  of  least  squares,  the  problem  becomes  that  of  investigating 
the  maximum  effect  of  random  errors  in  the  data  on  the  coefficients  of 
the  cubic.  This  is  done  in  a general  way  that  will  apply  to  a polynomial 
of  any  degree  and  a range  of  any  odd  number  of  stations. 

Curvature  changes  associated  wit a particular  distributions  of  errors 
are  found  to  cause  uhe  largest  errors  in  the  least  squares  coefficients 
and  hence  in  the  drag  values , 
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I.  INTRODUCTION 


1 i 

In  a recent  paper  Bergdolt  [l]^  gives  a description  of  the 
construction  and  instrumentation  of  a 12"  pressurized  range  for  small 
caliber  projectiles.  Some  preliminary  drag  measurements  (K^  vs  M)  are 

presented  and  it  is  noted  that  the  scatter  of  data,  particularly  in  the 
neighborhood  of  the  drag  curve  maximum,  is  unexpectedly  large . 

Analysis  has  been  made  by  Karpov  [2]  and  Charters  [3]  of  the  effect 

of  random  errors  in  time  and  distance  on  drag  data  obtained  in  firing 

ranges  and  in  combination  range -wind  tunnels  respectively.  To  simplify 

the  problem  both  authors  assume  that  a random  distance  error  e can  be 

s 

represented  as  an  equivalent  random  time  error  e^  by  the  equation 

e = v e.  where  v is  the  average  velocity  of  the  projectile.  The  treat - 
s t 

merit  of  Karpov  is  based  on  standard  statistical  methods  assuming  a 
least  squares  fitting  for  the  time  vs.  distance  data  and  a normal 
distribution  law  for  the  random  errors  at  each  station.  Charters  follows 
Karpov  in  this  respect  and  compares  the  two  types  of  facility  using  the 
Karpov  results.  Neither  author  attempts  to  determine  the  maximum  errors 
in  drag  caused  by  specific  distributions  of  error  among  the  stations  of 
the  range. 

In  this  report  we  find  upper  bounds  for  the  error  in  caused  by 

random  errors  in  a five  station  range  as  follows:  l)  by  examining  an 

idealized  problem  for  which  all  possible  distributions  can  be  easily 
enumerated  and  2)  by  establishing  bounds  for  any  possible  distribution 

of  errors  by  an  argument  based  on  the  results  for  the  idealized  problem. 

- • 

II.  STATEMENT  OF  THE  IDEALIZED  PROBLEM 


We  consider  here  the  case  where 

1)  the  number  n of  stations  is  small  e.  g.  n = 5, 

2)  errors  in  time  of  equal  size  g occur  at  each  station, 

3)  + or  - signs  are  equally  likely. 

For  convenience  in  calculation  we  consider  the  stations  equally 
spaced  and  the  origin  z «=  0 placed  at  the  central  station,  n is  always 
odd. 


III.  THEORY 

The  underlying  reason  for  considering  the  problem  stated  above  is 
the  intuitive  feeling  that  a certain  few  distributions  of  random  errors 
have  the  greatest  effect  on  the  least  squares  fit  of  time  t versus 
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Numbers  in  brackets  refer  to  the  bibliography  at  the  end  of  this  report 
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where  Z is  the  U x U matrix  above  and  a,  T the  corresponding  vectors. 
Since  the  Z matrix  depends  only  on  the  coordinates , which  remain  fixed, 
the  effects  of  a distribution  of  errors  d^  will  be  felt  only  in  the 

a and  T vectors.  For  true  time  i.  e,  distribution  c)  we  write 


(5)  Z a = T, 

and  subtracting  (5>)  from  (U)  yields 

(6)  Z (a  - a)  « T - T, 

Vfe  evaluate  T - T from  (1)  and  (3)  obtaining 

/ z° 


(7)  Z (o  - a)  - € 


1 

J '.j 


£ 


'V 


C. 

Z . 


I 


\ *> , 

From  this  we  see  at  once  that  the  numerical  error  A a in  any  coefficient 
is  proportional  to  € and,  once  the  Z matrix  is  determined,  dependent  only 
on  the  distribution  of  errors 

A further  result  may  be  obtained  more  heuristically  as  follows. 


For  a project  ire  or  u&oo  m and  uictftguBr  d Tired  in  atmosphere  of 
density  p 

(8)  K = 2 m a2 

a*  a., 

p d 1 

as  given  on  p,  8 of  [ 2^.  Comparing  rounds  fired  at  the  same  velocity 
vQ  = 1 but  different  densities  we  see  that 


( 9 ) a2  constant 

P 


must  hold  since  KQ  values  are  independent  of  density  - ignoring  for  the 
moment  small  Reynold’s  number  or  "scale"  effects. 

Thus  with  Z and  d^  given, 

*w 

do)  A a2  _ C (dk)  J_ 

a,  ~ 
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Thus  percentage  error  p?  in  a?  is  1)  proportional  to  and  2) 

c.  H 

dependent  on  a function  of  d,  , The  first  conclusion  agrees  with  Karpov's 
expression  for  p2  in[2j  middle  of  p.  Q.  The  latter  conclusion  is 
elaborated  below. 


IV . COMPUTATIONS 

To  test  the  dependence  of  errors  in  a-^  & a^  on  the  error  distribution 

dj.,  two  cubic  curves  were  chosen  as  representative  of  1 and  5 atmosphere 

firings  respectively.  The  coefficients  selected  correspond  t^  those  of 
least  squares  cubic  fitted  to  data  of  previous  experimental  drag  firings, 
modified  only  by  rounding  up  to  a convenient  number  of  figures, + Time 
values  were  computed  for  each  polynomial  at  z values  of  - 10',  - 5',  and 
0! . The  five(z,t)  points  are  then  called  the  true  values  and  correspond 
to  the  error  distribution  00000.  As  stated  in  § the  choice  of 
symmetrical,  equally  spaced  z's  ’.was  made  to  simplify  the  solution  of  the 

r •>  2n  1 

least  squares  normal  equations,  since  with  this  choice  = 0 

for  any  integral  n >0.  The  least  squares  cubic  for  the  5 atmosphere 
case  was  determined  for  each  of  the  16  error  distributions  by  a Gaussian 
elimination  scheme.  A change  in  sign  of  s is  equivalent  to  reversing  the 
sign  of  the  error  distribution  which  yields  a percentage  error  in 
2a 

2 of  the  same  magnitude  but  of  opposite  sign,  As  a check  on  the 
al 

computations  this  was  verified  for  k such  distributions. 

The  components  of  th«2  error  vectors,  (a  - a),  were  computed  by 
solving  ( f) . Since  (a  - a)  is  dependent  only  on  the  z's  and  e , it 
applies  to  both  the  1 and  c.  atmosphere  cases  in  the  determination  of 
the  percentage  error  in  2a2/a^ , 

V.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

The  main  results  of  th°  computations  are  summarized  in  T~ble  I 
which  gives  the  16  error  distributions,  their  corresponding  (a  - a) ' s 
and  percentage  error  in  K~,  for  the  1 and  5 atmosphere  cases.  Reversing 
every  sign  in  Table  I will  give  the  data  for  the  other  half  of  our 
(m  .)  matrix. 


The  choice  of  e = 1 is  justified  as  follows:  Our  e equals  Karpov's 

et,  see  C2j  page  8,  where 


(11)  e 


2 2 

t = et  et 

t t u2 


and  e+  , and  e.  are  time  and  equivalent  distance  errors  respectively. 
”1  2 


3 


r 


1 
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The  least  count  of  our  counters  is  0,625  p sec.  and  the  estimated  error 
in  any  z determination  is  .001  ft.  v.  [i] . If  we  let  e.  * 0.625  u sec, 

with  e^  = £ =1  then  e^  = .76  which  corresponds  to  a distance  error 

of  .0025  ft.  at  a velocity  of  3000  feet  per  second.  It  seems  unlikely 
that  our  distance  error  is  this  large;  therefore  the  £ = 1 value  might 
reasonably  be  somewhat  diminished.  For  convenience  it  was  not. 

In  the  5 atmosphere  case,  the  maximum  percentage  error  in  Kn, 

.91#  occurs  for  distribution  (l5)  + -+-  - +.  Hence  all  nercentage 
errors  will  be  contained  in  the  band  - .91#  about  the  true  curve.  If 
all  distributions  are  considered  equally  likely,  this  extreme  value 
occurs  in  l/l6  of  the  cases.  Two  other  distributions  causing  curvature 
change,  (?)  * + - - + and  (13/  +--++,  have  percentage  errors  of 
»7h#  and  .63#  respectively.  All  others  (i.  e.  13  remaining  distributions) 
have  absolute  percentage  errors  of  .1*6#  or  less . 

Similarly,  the  error  band  for  the  1 atmosphere  case  is  found  to  be 
- i**60#  or  five  times  the  width  of  the  5 atmosphere  case.  13  distributions 
out  of  16  lave  absolute  percentage  errors  of  2.20#  or  less. 

If  the  scatter  observed  in  the  experimental  data  is  due  to  the 
assumed  random  errors  of  this  idealized  problem,  it  will  be  bounded 
by  the  above  limits.  If,  however,  points  fall  outside  these  limits, 
then  1)  the  estimated  £ is  too  small,  2)  errors  other  than  random 
timing  errors  contribute  to  the  scatter  or  3)  the  analysis  of  the 
idealized  problem  does  not  apply.  We  show  later  that  the  bounds  on 
error  in  for  any  distributions  such  that  jin^  4 1 differ  only 

slightly,  if  at  all,  from  those  already  obtained.  Thus  the  idealized 
problem  furnishes  reasonable  bounds. 

v 

The  scatter  observed  in  the  vs.  M curves  in  Bergdolt's  report 

[1]  for  1 and  5 atmospheres  is  bounded  by  the  limits  mentioned  earlier, 
if  a few  extreme  points  are  disregarded.  There  are  several  valid 
reasons  for  disregarding  this  small  number  of  extreme  points.  At  the 
time  the  rounds  represented  in  Figures  9>  10  and  13  of  £l]  were  fired, 
there  was  no  control  of  the  humidity  of  the  air  in  the  range  and  no 
system  of  checks  in  the  operation  of  the  counters.  Later,  it  was  also 
discovered  that  the  calibration  of  the  pressure  gauge  wa*  in  error. 

Any  of  these  three  types  of  error  might  cause  points  to  fall  outside 
the  limits  stated  above j however  subsequent  experience  with  the  counters 
indicates  them  as  the  likely  source. 

Since  these  three  difficulties  have  recently  been  minimized  by  the 
use  of  dry  air,  auxiliary  check  counters  and  recalibrated  pressure  gauges, 
we  shall  expect  more  recent  data  to  fall  within  the  bounds  specified. 

In  fact,  returning  to  Table  I,  it  should  be  noted  that  13/16' s <j)f  the 
distributions  have  percentage  errors  falling  wit^an  the  bounds  - 1/2 
maximum  percentage  error  i.  e.  for  1 atmosphere  - 2.20#  and  for  5 
atmosphere  - .1*6#.  Hence  we  may  conjecture  that  in  the  practical  case, 

80#  of  the  points  will  cluster  within  the  central  half  of  the  maximum 
band  width. 


i 
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The  combined  counting  and  distance  errors  nay 

V 


6 1. 


in  actuality  take  on 
only  lie  in  the  range 


We  prove  in  Appendix  1 that  all  distributions  with  m^j's  which  fall 

within  this  range  have  errors  in  bounded  by  limits  only  slightly 

larger  than  those  just  obtained  for  m^  = - 1.  It  seems  likely  that 

these  latter  are  also  true  bounds  for  any  distribution  such  that 
j rn^_.|  sj  1.  'While  a proof  of  this  conjecture  has  not  been  obtained,  no 

exception  has  yet  been  found.  Its  truth  or  falsity  is  not  of  much 
importance  at  present  since  the  limits  actually  obtained  in  Appendix  1 
only  differ  by  a small  positive  quantity  from  those  for  m = - 1 , 

(See  Table  II) 


Thus  we  may  conclude  that  the  scatter  band-width  indicated  by 
our  idealized  calculation  for  the  32  distributions  obtained  above 
gives  an  accurate  estimate  of  the  band-width  to  be  expected  in  practical 
experience  with  the  five  station  range.  It  is  to  be  understood  that 
suitable  adjustments  in  the  assumed  £ must  be  made  as  changes  occur  in 
the  fundamental  accuracy  with  which  time  and  distance  can  be  measured. 


F.  D.  BENNETT 
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APPENDIX  1 


L 


The  computations  discussed  earlier  are  limited  to  an  £ of  1 and 
of  either  - ] 

-1  S ^ + 1 , 


m,  . of  either  - 1.  Now  we  shall  examine  other  possibilities  for  which 
KJ 


From  § HI  we  have 


\ - / — 


2m 


P d 


By  logarithmic  differentiation  we  get 


(A-l) 


A a. 


2 dl 

It  is  easily  shown  that  and  a^  are  both  positive  quantities 
Therefore 


(A -2) 


AKr 


I A a2 1 . I 


A 


Now 


(A-3) 


Aai 

A a0 
Furthermore 
al 


/ \ i \ 

vA-u; 


(A  -5) 


(A-5a) 


Thus, 

Aiq 


D 


I Aal  I 

Aa2  | 


max 


max  . 


" al 

+ Aa^ 

31  ” 

H 

n) 

<1 

max 

“ *2 

♦ A 

> 

i 

i 00 

Iw 

A a2 

max 

Aki 

■v 


^ 

Aa2j 

max 

Aal 

max  or 

max  a?  - 

* A2  « 

Aa2 

. A 

1 

a,  - Aa, 

max  1 I 1 

• 

max 

max 
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In  order  to  find  the  bounds  on 


it  is  now  necessary  to 


find  values  for  ! A a,  I and  j A a,, 

Fror»  § m we  nave 

2(a  - a)  - « / \ • 

. ri  \ 


AKD 

V 

MX 


(7) 


y m 2 
^"kd  zd 


\smkd 


5d 


We  solve  for  the  individual  Aa^'s  and  obtain  the  following 


expressions 
(A -6)  A a, 

(A-7)  A a. 


g |-ML2  Jo  + **22  J 1 

"deFT  [-Mjg  2v3 


and 


6 

det  t 


f +M 

V* 


X „ - m X „ 


”l3  ^3 

S3  ^2  - % 


1 

► 

3 


where  X.,  » Xi*..,  z^j  and  M.  ^ is  the  signed  minor  or  cofactor  of  tbs 


.th 


H.J 


‘id 


i .i — element  of  the  Z matrix. 


For  our  particular  Z matrix,  all  M.  .*s  for  which  (i  + j)  = 2n  + 1 


are  equal  to  zero.  M «s  for  which  (i  ♦ j)  * 2n  are  constants  which  we 


can  evaluate. 

V 

Hence  we  may  simplify 

(A-6)  and 

(A-7)  as 

follows 

(A -6a) 

A 

6 Tiu.  x. 

• fc  j CC  JL 

det  Z L 

r~ 

M.^ 

4<- 

M 

U> 

J L — 

(A -7a) 

A ag 

- 6 Pl3  XQ 

det  2 

+ 

’b) 

To 

find 

Aai  mxand  |Aa2 

max 

we  examine  the 

behavior  of  the 

bracketed  functions  as  follows. 
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For  Aa-j: 


F‘22  ^1  * 'V  ^*3  " ‘*22  2j  * m1j2  ^“kj  z3 


Smk3  [fl2 


b2j  * * 


ti 


( ii  - Kgg  2j  * "1,2  2,>- 


Analysis  of  . for  j - - 2,  - 1 and  0,  summarized  in  Table  III, 


shows  that 

O 

the  maximum 

l*\i 

+3 

-1,  tl,  *1, 

-1  and  -1, 

♦I.  -1 

, -l. 

must  yield 

a smaller  | 

2"kj 

*jl 

Siiailarly  1 

'or  Aa^i 

Vi13  ^ 0 + 

H33 

•i3 

Zm 

" 

[“ 13 

0 
z . 
0 

= S"kj 

0j  '• 

% ' 

-L  _ A m,  . z . 

33  icj  j 


Distributions  giving  maximum  |Z}m^  j are  +1,  -1,  -1,  -1,  +1  and 

-1,  +1,  +1,  +1,  -1,  (see  Table  III).  This  analysis  is  verified  by  the 

computations  summarized  in  Table  I.  It  may  be  noted  that.  lAa.l 

i d i max 

and  | A a,  I do  not  occur  simultaneously.  Hence,  we  conjecture 
1 1 nax  ^ ■ + 

that  for  a given  e value,  11  is  bounded  by  - a . Ir.  practice 

KD 

a K 

the  contribution  of  the  a.,  term  to  D is  small;  and  may  be 

KD 

neglected  so  that  the  preceding  expression  gives  a good  estimate  of  the 
bound  * 


For  actual  computation,  any  common  numerical  factors  in  jr  . were  factored 
out  and  combined  with  e , ^ 

det  Z 


w- 
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TABLE  I 


Aao 

Aal 

Aa2 

Aa3 

(1)+++++ 

1.00 

0 

0 

0 

0 

0 

(2 )++++- 

1.17 

.03 

-.011 

-.0013 

-.1*6 

-2 .20 

(3)+++-+ 

.31 

-.27 

.006 

.0027 

.28 

1.25 

(1*)+++ — 

.1*9 

-.23 

-.006 

.0013 

-.18 

-1.16 

(5 )++-++ 

.03 

0 

.011 

0 

.1*6 

2.20 

(6)++-+- 

.20 

.03 

0 

-.0013 

-.01 

-.01 

(7)++—+ 

-.66 

-.27 

.017 

.002? 

.71* 

3.1*5 

(8)++ — 

-.1*9 

-.23 

.006 

.0013 

a r* 

1 .?!;  i 

(9)+-+++ 

.31 

.27 

.006 

-.0027 

.18 

1.11* 

(10)+-++- 

•k9 

.30 

-.006 

-.001*0 

-.29 

-1.26 

(11)+-+-+ 

-.37 

0 

.011 

0 

.1*6 

2.20 

(12 '»+-+— 

-.20 

.03 

0 

-.0013 

-.01 

-.01 

(13)4*— 

-.66 

.27 

.017 

-.0027 

.63 

3.3U 

(1U)+— +- 

-.ii9 

.30 

.006 

-.001*0 

.17 

1.11* 

(15)+ — ♦ 

-1,3<! 

o 

.023 

0 

.91 

A*.6o 

(16)+ 

-1.17 

.03 

.011 

-.0013 

.1*5 

2.19 

5 Atmosphere  Time  Curve 


t « 5000  + 500  z + 2.5z2  + .02z3 
1 Atmosphere  True  Curve 


t = 3800  + 520  z + .5z2  + .01  z3 
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